NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



• ■& <C-~ if* 




1 RESEARCH ON OUTPUT FEEDBACK CONTROL 


INTERIM REPORT 


January 1985 - September 1985 
November 1985 


Research Supported by NASA-Langley Research Center 
NASA Grant No. NAG-1-243 


Principal Investigator: Dr. Anthony J. Calise 

Research Assistant: Mr. Friedrich S. Kramer 

NASA Grant Mont tor: Dr. Douglas B. Price 


MECHANICAL ENGINEERING AND MECHANICS 
DREXEL UNIVERSITY 
PHILADELPHIA, PA 19104 


{NASA-CP- 174253) 
FEEDBACK CCHTBOL 
Sep. 1985 {Drexel 


RESEARCH OH OUTPUT 
Interic Report, Jan. ~ 
Unit.) 20 F HC A02/BF 

CSCL 


A01 

12A 


1186-13028 


Unclas 

047 86 





G3/64 


SUMMARY 


The research for the current year has dealt with the synthesis of fixed 
order dynamic compensators for multivariable linear time Invariant systems by 
minimizing a linear quadratic performance cost functional. In particular, 
attention has been given to robustness Issues In terms of frequency domain 
specifications, namely low frequency performance, high frequency robustness to 
unmodeled dynamics and parasltlcs, and cross-over frequency sensitivity. 

In designing fixed order compensators, an output feedback formulation has 
been adopted by suitably augmenting the system description to Include the 
compensator states. However, the minimization of the performance Index over 
the range of possible compensator descriptions was Impeded due to the 
nonualqueness of the compensator transfer function. A controller canonical 
form of the compensator was chosen to reduce the number of free parameters to 
Its minimal number In the optimization. In the MIMO case, the controller form 
requires a prespecified set of ascending controllability Indices. This 
constraint on the compensator structure Is rather J .r..i^<'uous In relation to the 
increase In convergence rate of the optimization. Moreover, the conti oiler 
form Is easily relatable to a unique controller transfer function description. 
This structure of the compensator does not require penalizing the compensator 
states for a nonzero or coupled solution, a problem that occurs when following 
a standard output feedback synthesis formulation. 

Unlike full state LQR designs, output feedback designs do not have 
guaranteed gain and phase margin properties. For that matter, neither does 
an LQG design. Recent robustness recovery techniques try to asymptotically 
approximate the full state LQR design by Incorporating fictitious noise in the 
performance Index. A cheap control formulation results with high frequency 
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dynamics appearing in the compensator. The approach taken in this research is 
to extend the concept of frequency-dependent cost functionals to the output 
feedback problem. Lead-lag performance weightings have been used to 
manipulate the performance index and provide a method to expand the 
compensator to enhance the systems closed loop robustness. Results on SISO 
systems confirm our hypothesis that output feedback dynamic compensators can 
be designed to meet phase margin requirements of the system without 
sacrificing control bandwidth constraints. Design examples for MIMO systems 
are currently underway, and similar results are expected. 

This report summarizes the results to date on the SISO case. The 
contents of this report will be submitted for pr^ ntation at the next AIAA 
Conference on Guidance, Navigation and Control. 
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1. Introduction 


It Is well known that Linear Quadratic Regulator (LQR) synthesis methods 
have guaranteed stability margins. Unfortunately, these properties hold 
only in case of full state feedback. Observer based compensator design 
techniques exist to estimate unavailable plant states, and make the LQR 
design viable. However, this combination of state estimation and regulation 
may result in a compensator with poor stability margins, even though the 
separate designs are robust. It has been shown that the robustness properties 
of the LQR design for nonminimum phase plants can be recovered via an 
asymptotic method called Loop Transmission Recovery (LTR) [1]. The LTR 
method relies on a cheap control formulation with a subset of the compensator 
dynamics becoming fast. In this post analog era, finite word length, sampling 
rate or time delay may impede the use of high gain controllers when 
implemented digitally. Several other shortcomings of the LTR method that are 
related to plant inversion are discussed in f 2 1 . Aside from the robustness 
issue, the order of the compensator when designed for large scale systems may 
prove unwarranted. 

Optimal output feedback design of fixed order compensators introduced in 
the early seventies [3] has received limited attention throughout the years. 
Part of the difficulty rests with the inability to characterize the stability 
margin properties. Unlike the algebraic Ricatti equation that arises in the 
case of full state feedback, the necessary conditions that result from the 
optimal output feedback problem are not conducive to analysis in the frequency 
domain. The other difficulty in synthesizing optimal dynamic compensators is 
that the standard approach of adjoining the compensator dynamics to the plant 
dynamics (and reformulating as a static gain output feedback design) results 
In an over parameterized formulation. This is a direct consequence of the 


fact that the compensator lacks a predefined structure, which Invariably 
results In difficulties with convergence to an optimal solution. This Is 
avoided here by representing the compensator in controller canonical form. 

This paper presents an LQ optimal design approach for designing fixed 
order compensators that satisfy frequency domain performance and robustness 
specifications. Only single-input single-output plants are considered; 
however, the approach can be extended to the multivariable case. Implicit in 
the system description is the presence of noise in the measurement, and to 
avoid control activity gent v ^d by this noise direct output feedthrough in 
the compensator has been eliminated. 

The following section discusses the optimal outpu': feedback formulation 
in controller canonical form. Next, it is shown how the concept of 
frequency-dependent cost functionals can easily be handled in this canonical 
setting. This is useful in improving the stability margins of the 
compensator. Finally, examples are given to elucidate the design technique. 


2. Fixed-Order Dynamic Compensation 


Consider the linear 

time invariant 

single-input single-output 

system 

• 

x - Ax + bu 

x(0)«x o 

x t R n 

( la) 

y ■ c^X 


y e R 1 

(lb) 

The system (1) is controlled by the fixed-order dynamic compensator 

without 

direct feedthrough of the output 



u ■ -h^z 


u e R 1 
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z * -Pz - ny 

z(0)*0 

z £ R n c 
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where 0<n c <n Is chosen a 

priori subject 

to the requirement that there exist 

stabilizing compensator 

for the closed 

loop system defined by 
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The dynamic compensator (2) is defined In a controller canonical form: 


u = -h t z 
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(4a) 
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h fc - [hi . . . hn ] p* - [pi . . . Pq ] (5) 
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This canonical form reduces the number of free parameters from n c ^ + 2n c to 


■*» 


j 

I 

i 

2n c . Furthermore, the companion form of P ■ P°-np t provides structure to the 
compensator and Is a natural form for eliminating direct feedthrough of the 
t output. This point will be further explained below. 

In the spirit of an output feedback formulation, augment the plant states 

l ^ /V 

' to Include the compensator states, that is, let x* ■ [x c , z t ], y c - [z 1 ] and 

u* ■ [u,u c ]. The augmented system Is then given by 


x * Ax + Bu 
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Note that the controller form reduces the number of compensator control 
elements (u c ) from to 1. The fixed order dynamic compensation problem has 
now been converted to a static output feedback problem where the compensator 
gains (h,p) are given by G. Notice that this formulation alleviates the need 
for zeroing those elements associated with the static feedthrough which would 
arise If the standard output feedback formulation is used. 

The compensator transfer function 


( 8 ) 


K( 8 ) - h^slj, +P) -1 n - h(s)/p(s) 

n -1 n -1 n 

- (hi + ... + h„ s c ) / ( P i + ... + Pn s c + s c ) 
c c 

Art 

directly correlates with the elements of G in (7). Note that the controller 
form is easily relatable to a transfer function description, and as such any 
parametric optimization of individual compensator gains provides a direct 
transfer function counterpart. In particular by zeroing the ascending 
coefficients Pi,P2».** of the monic polynomial p(s) defined in (8), type 
1,2,... feedback systems are formed. Likewise the degree of compensator 
rolloff can be controlled by zeroing the appropriate elements in h 1 . 

The resulting composite system can now be ast as a parametric optimi- 
zation problem. The free parameters are chosen to minimize the 

quadratic performance index 

OD ~ <«•*» 

J(h,j>) >= E{ / [x c Qx + pu^ + p c u c 2] dt} (9) 

o 

where E denotes expectation over the distribution on x (0), Q * D l D such that 

(A,D) is detectable, and p and p c are positive scalars. The input u c « -g^z-y 

drives the compensator, and increasing or decreasing the penalty on u c 

correlates directly with decreasing or increasing the compensator bandwidth. 

The solution to the LQ optimal output feedback problem depends on the 

initial distribution of the states. Extending the standard LQ c tlmal output 

feedback design to Include a compensator requires some knowledge of the 

compensator states' initial distribution. One approach is to assume that the 

compensator states are known and at reBt, and that E{[x 0 ,z 0 ] [x^Zq] 1 } “ block 

diag [X^jjZq], Zq - 0q . In practice however, the initial distribution of the 

c 

compensator states may be needed for the existence of (CLC^)~^ in the calcu- 
lation of G. This fact 1 8 stated in the following thereon, which is proven in 
the appendix. 
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Theorem 


A sufficient condition for the existence of (CLC*)" 1 with (P ,rt > 

In controller form is that Z„ 3 nn t . 

3. Frequency Shaped Conpensation 

Frequency-dependent cost functionals have been recently Introduced as a 
way of embedding classical design concepts within the context of LQ optimal 
control [A, 5]. These frequency-dependent cost functionals have been developed 
for the design of a full state controller and Its dual, the state estimator. 
When applied to the design of observer based compensators, the resulting 
controller Is dimensionally larger than the system model due to the 
introduction of the frequency shaping terms [4], In this section we formulate 
the use of frequency-dependent cost functionals In the context of LQ optimal 
output feedback. 

Modifying the performance Index In (9) to Include Input and output 
frequency shaping, the cost functional becomes 

J(£,h) - E{ l/2it / [x*( ja,)Qx(Jo,) + |Q y (jfc)y( ju-) | 2 + p |R( j<Ju( jo,) | 2 + 

PcMJ<“>l 2 l<k} (10) 

where * denotes complex conjugate transpose, Q, p, and p c are as before, and 
Qy(j(L> and R(ja,) are frequency shaping terms with the particular lead-lag 
form: 

Qy ( j u, ) ■ (yjo, + l)/(0jou + 1) (11a) 

R(jcJ - (Bjo. + l)/(ajo* + 1) (11b) 

The weighting functions Q y (j (1 ,) and R(jt) were chosen for their effectiveness 
In shaping the stability margins In classical designs; however, other rational 
functions of frequency are possible [4]. Realistically one would select B > a 
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> 0 to penalize higher frequency control activity and 6 > y > 0 to weight low 
frequency sensor Information. More Importantly, the frequency shaping 
parameters are varied to adjust the loop phase and gain margins while still 
optimizing the compensator. 

The cost functional In (10) can be expressed in the time domain as 


J(£,h) “ £ (xtQx + yj 2 + puj 2 + p c u c 2 }dt (12) 

with the system (6,7) augmented to Include the state realization of the 
frequency shaping terms (11): 


• . 

wj - wj/e y W[ e (13a) 

2 

yi * wi(e-nr)/Q + y/0 yi e Rl (13b) 

W 2 * W 2 /B + uj wj c R* uj e R* (14a) 

2 

u * W2(B~a)/B + u^a/B (14b) 

The dynamic compensator Is expanded to include the states wj and W 2 : 

U 1 “ " 811 w l ” 812 W 2 ~ h tz (15a) 

z ■ P°z + nu c ( 15b) 

u c “ - 821 w l * 822 w 2 ~ £ tz “ yi 05c) 


Design via Optimal Output Feedback 

The optimization problem can be reformulated using the augmented system 

^ a, e. 

description (6) where x fc ■ [x 1 , wj , W 2 , z fc ], y l - [wj, W 2 , z t l, u c “ [uj, u c l 
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The frequency shaped fixed order compensator design is now in terms of a 
static gain output feedback design, and the solution G does not require any 
zeroing of elements. However, as ;:ill be shown in the next section, zeroing 
the gn and g22 terms in G gives a classical interpretation of the effect 
of frequency shaping in the case of output feedback. 


Transfer Function Description 

The transfer function description of the fixed order dynamic compensator 
with lead-lag frequency shaping (13-15) is given by 

u(s)/y(s) * K(s) ” (as + l)K' c (s)/(es +1) (17) 

where 


K' c (s) ■ 1 (YS + l + g 21 Y)Kc(s) " gliy]/( (f5s + l + f?l2fl)”822B K c( 8 )1 (18a) 


K c (s) ■ h t (eI+P)“^ii 


(18b) 


By zeroing gn and § 22 * the frequency shaped compensator, K(s), 
simplifies to 

K(s) - F u (s)K c (s)F y (s) (19) 

where 

F u (s) " (aS+l)/( 6 s+l+gi 2 B) (20b) 

Fy(s) ■ (YS+l+g 2 iY)/( 6 s+l ) (20b) 

which Is lust the nominal compensator (with Its Internal parameters 
reoptimized) cascaded fore and aft with the lead-lag type first order filters, 
F u (s) and F y (s), respectively. 

Both the Input and output frequency shaping cost functions (10) were 
Introduced to point out the benefits of each. Notice from (17) that by 
shaping the output of the system with a lead-lag cost function Qy(!k) In 
(11a), an pole Is placed In the compensator at -1/e. The location of the zero 
In the compensator is Influenced by the gains g 2 j and 811* Conversely, by 
shaping the Input of the system with a lead-lag cost function R( j a.) In (lib), 
a zero la placed In the compensator at ~l/a. The location of the pole In the 
compensator is Influenced by the gains gj 2 and g22* Hence, because of the 
commutative property of SISO systems, the only difference between shaping the 
output or Input of the system Is that either a pole or zero of the resulting 
compensator Is fixed by the weighting parameters 6 or a. 

4. Example 

Consider the simple 2nd order SISO system with state description: 


X 1 “ X2 

X 2 ■ -3xj - 4x2 + u 
y - 2xj + x 2 


( 21 ) 


and loop transfer function: 

G(s) - (s+2)/((s-rl)(s+3)l (22) 

The performance specifications call for 20 dB of gain for frequencies 
less than 1 rad/s and uiross-over frequency 10 rad/s. 

A first order compensator (n c »l) was designed to meet the performance 
specifications with associated cost penalties 

0 ■ block dlag[Q x ,10^1 Q x - dd t , • 4/ 5 ( /35 , 1 1 10~2 
p - 0.0035 p c - 25.0 (23) 

and 

[X 0 ,Z 0 ] - diag[ 1,1,1] (24) 

The example and plant weighting matrix, 0 X , were taken from [6], where an LTR 
design was Illustrated. The algorithm described In [7] was used to solve the 
optimal output feedback problem. This algorithm has guaranteed convergence 
properties. It also permits a constraint of the form y(G) « 0 to he 
approximated using a penalty function approach. This can be used to zero 
selected elements In G, or to place selected poles or zeros o' the compensator 
as will be illustrated below. The above performance weightings resulted In 
the following compensator design 

K(s) - 124.9/( s+4.86) (25) 

Next, a second order compensator was designed to exhibit type 1 behavior 
and have a zero at 1 rad/s. This design was carried out by zeroing the pi 
element In G of (7) and by constraining hj to equal h2- The cost penalties 
were the same as above with Q and Z 0 expanded to Include the extra compensator 
state 

Q - block dlag [Q x ,0,103] and IX 0 ,Z 0 ) - dlag 11,1,0,1] (2o) 

The performance Index cost Increased by 15Z and resulted in the compensator 
K(s) - 129.9(s+l)/s(s+5) (27) 
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The magnitude and phase responses of the loop transfer function G(s)K(s) for 
the first order compensator (25) and the type 1 second order compensator (27) 
are Illustrated on a Bode diagram in figure 1. 

Lead-Lag Augmentation 

The preceeding first order compensator design (25) meets the performance 
specification defined earlier in this example and compares favorably with the 
full state feedback controller designed in (6]. The first order compensator 
results in 35 degrees of phase margin; while, the full state design gives 85 
degrees. Manipulating the input and output performance weightings to improve 
the phase margin of the first order compensator resulted in designs which 
strayed from the original specifications with no reel improvement in the phase 
margin. The performance index was next augmented to include a lead-lag term 
at the output 

Q y (j(J - (0.2ju. -*-l)/( .02ju, + 1) (28) 

which from classical compensation techniques should add about 50 degrees of 
phase margin near 10 rad/s. rhe compensator was then redesigned with the 
associated cost penalties 

Q « block dlaglQ x , 0,10^] , p ■ 0.30 p c ■ 40.0 
and (29) 

lXo,Z 0 ] - diag [1,1,1,11 

Frequency shaping resulted in the compensator transfer function 

K(s) - 592.4(s+4.28)/(s+2.05)(s+50) (30) 

The phase margin was increased to 75 degrees. Bode diagrams of the loop 
transfer function G(s)K(s) for the first order compensator (25) and the second 
order compensator with frequency shaping (30) are compared in Figure 2. 
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5. Conclusions 


A design methodology has been introdi ;ed for the LQ synthesis of fixed 
order dynamic compensators without output feedthrough for SISO systems. A 
controller canonical structure was imposed on the compensator description 
which minimized the number of free parameters and correlated with a unique 
transfer function realization. A frequency-dependent cost functional provided 
the means to enhance the robustness of the fixed order compensator design to 
improve the stability margins of the closed loop system. In particular, 
first order lead-lag type cost functionals were introduced to shape the input 
and output of the system and were shown to have a classical lead-lag 
counterpart when used with the dynamic compensator. 
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Appendix 
Proof of Theorem 


The necessary conditions satisfying the LQ optimal output feedback 
problem are given by 

A fc K + KA + Q + C^tRGC - 0 (Al) 

AL + LA 1 + XX 0 « 0 (A2) 

RGtCLCt) * BtRLC 1 (A3) 

a 

where A * is asymptotically stable, the quadruple (A,%,£,G) is defined 

in (6) and (7), Q > 0, R * diag [p,p c ] > 0 and XX 0 ■ block diag [X 0 ,Z 0 ]> 0. 

<"• /s. 

In order to uniquely determine G in (A3), it is necessary that (CLC^) -1 exist. 
For the form of C in (7), this inverse can be reduced by first partitioning L 
into 


( A4) 



L 11 

l 12 

L * 

t 



_ L 12 

L22_ 

and expanding CLC t . 

Thus 


CLC c - L 2 2 

and the existence of (CLC) - * is solely dependent upon the existence of 
(L22)” 1 . 

Expanding equation (A2) results in 


(A5) 


PL 2 2 + L 22 P^ + Iq ? cn ^ + nc^hj 2 + Zq * 0 


(A6) 


- 15 - 


Referring to [8] and the references therein, the following remarks hold. 

Equation (A6) can be written as 

+ + 

(P + + L22(P + **£^^ , 12^ , 22^ + Z 0 * 0 (A7) 

where + denotes a pseudoinverse. Since (P,n) is controllable, then so is 
+ 

(P + nc t Li2b22»“) • Let z e " on*. Furthermore, since L22 is at least 

+ + 
non-negative and (P + nc^L; 7L79 ,n) is controllable, then (P + nc t L;2f'22) is a 

stability matrix; and the solution L22 to (A7) is positive definite. 

+ 

We can broaden the class of Z 0 to that which makes (P + nc^L; 2L22 ,D) 
controllable, where DD fc * Zq 3 nn fc . 
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